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We investigate the concept of a standard map for the interaction of relativistic particles and 
electrostatic waves of arbitrary amplitudes, under the action of external magnetic fields. The map 
is adequate for physical settings where waves and particles interact impulsively, and allows for a 
series of analytical result to be exactly obtained. Unlike the traditional form of the standard map, 
the present map is nonlinear in the wave amplitude and displays a series of peculiar properties. 
Among these properties we discuss the relation involving fixed points of the maps and accelerator 
regimes. 



I. INTRODUCTION 

Wave-particle interaction has always attracted a great 
deal of attention as an efficient way to particle accel- 
eration and particle heating. A wide range of applica- 
tions indeed flourishes, from heating and current drives 
in fusion devices [H, 0] , to more alternative concepts like 
plasma based particle accelerators [3| and non-neutral 
beams Q. 

Wave-particle interaction is basically a nonlinear pro- 
cess [a, Q- This means that in the usual scenarios of 
the interaction one can expect to see regular and chaotic 
patterns interwined in the appropriate phase-spaces Q. 
Regular regions are useful for coherent acceleration while 
chaotic regions are adequate for particle heating. The 
prevalence of one or another type of region is a direct 
result of the strength of the perturbations impinging on 
the the particle motion. 

As shown over the years, investigation of nonlinear 
wave-particle interaction can be largely aided if area pre- 
serving Hamiltonian maps can be constructed for the sys- 
tem under study. Area preserving maps provide such a 
powerful tool in the study of the nonlinear dynamics of 
Hamiltonian systems because they yield a series of exact 
analytical results to work with [3, Q • 

Among the large variety of area preserving maps, the 
most prominent is perhaps the standard canonical map 
in one-degree-of- freedom and its close variants p, \W^ . 
The standard map describes a large number of systems, 
among which one can find the case of wave-particle non- 
linear interaction where particles move under the action 
of electrostatic fields [a, Q . Under these circumstances 
the map clearly shows how appropriate resonant condi- 
tions give rise to efficient acceleration mechanisms. 

Maps can be constructed from the underlying Hamil- 
tonians as one makes the assumption that the system 
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is periodically perturbed by impulsive kicks, followed 
by integration of the resulting canonical equations over 
one perturbing period. Impulsive approximations ade- 
quately represent pulsed perturbations, resulting from 
broad band spectra of count er prop agating waves driven 
by nonlinear wave coupling [ill . |l2| . 

In the standard case, while one of the canonical vari- 
ables appears only in the interaction term of the origi- 
nal Hamiltonian, the other appears only in the free part. 
This peculiar feature allows to show that one of the vari- 
ables is constant between the kicks and the other does 
not vary across a kick, which is the key factor to simplify 
the final form of the theory. As a further result, the map 
can be shown to be linear in the wave amplitude. If one 
now takes another case of interest where particles move 
under the action of an external magnetic field, the canon- 
ical variables become entangled by the field and can no 
longer be easily segregated into the free and interaction 
terms of the corresponding Hamiltonian. However, as 
will be shown, a set of convenient variable changes allow 
for the construction of an exact map. 

We shall focus on electrostatic perturbing modes prop- 
agating perpendicularly to the external magnetic field. 
These modes are frequently present and can be respon- 
sible for a large amount of particle energization. In 
open boundary systems, traveling electrostatic waves 
have been shown to produce efficient coherent acceler- 
ation [l|, [ij, [ij] as well as incoherent heating [T^. In 
closed systems, stationary modes formed with counter 
propagating waves are the ones to be considered. Closed 
systems encompass the cases of magnetically confined 
plasma columns and beams, and support waves with 
frequencies in the vicinity of the cyclotron frequency 
[la . Il7| . Resonant cyclotronic interaction is also possible 
between lower hybrid modes and fast electrons when the 
electronic cyclotron frequency decreases due to relativis- 
tic effects [18] . Large resonant islands [19] will be shown 
to form due to relativistic nonlinearities of the transverse 
particle motion with respect to the magnetic field, a fea- 
ture of relevance in accelerator regimes as discussed later 



The central interest of the present paper is then to 
construct the exact map for the dynamics of magnetized 
particles under the action of perpendicular electrostatic 
waves. The resulting map could be seen as the mag- 
netized counterpart of the classical standard map. Al- 
though related, we shall see that both maps display very 
different structures: nonlinearities in the dependence on 
the wave amplitude affects positioning of fixed points 
and the associated accelerating regimes of the magne- 
tized case. Limits on the acceleration efficiency of the 
magnetized case are discussed as well. 

The paper is organized as follows: §2 is devoted to the 
construction of the map, §3 to the analytical and numer- 
ical investigation, and in §4 we draw our conclusions. 



II. FROM THE MODEL TO THE NONLINEAR 

MAP 

Consider a particle with charge q, mass m and per- 
pendicular canonical momentum pj^, moving under the 
combined action of a uniform magnetic field B = 
BqZ and a stationary electrostatic wave of wavevec- 
tor fc, period T and amplitude A, lying along the x 
axis. Working with a set of dimensionless quantities, 
p±/mc — > p^, qBox/mc — >■ x, {qBQ/m)(t,T) — >■ (t,T), 
{\ / mc?){qBQ / m) A — > A, the dimensionless Hamiltonian 
H — > H /rac^ governing the particle motion can be writ- 
ten as 



H= ^1 + pI + x"^ + Aco&{kx)^5 {t - nT) , (f) 

n 

where with no loss of generality we consider the canonical 
momentum py = 0. We point out that even though the y- 
component of the canonical momentum is conserved and 
taken to be zero, dy/dt is not zero and the motion is not 
one-dimensional, c is the speed of light and, as stated ear- 
lier, for our purposes we focus on pulsed systems whose 
action is represented by the periodic collection of delta 
functions. 

At this point one must adopt a strategy to integrate the 
dynamics generated by the Hamiltonian ([I}. To this end. 



we note that particles are subjected to kicks whenever 
the time is a multiple of the period T. Between consecu- 
tive kicks, however, the one-degree-of-freedom dynamics 
is time independent and integrable, thus representable 
in terms of action-angle variables Px = V2Icos9, x — 
\/2Ism9 with constant action /. For future purposes, 
one can promptly find the form of the Hamiltonian writ- 
ten in terms of action-angle coordinates. It reads 

H = VTT2I + A cos{kV2I sine) J2Ht~nT) . (2) 



Across the kick, both / and 9 undergo abrupt varia- 
tions since both are present in the wave-particle coupling 
term. However, in the original variables x, Px only px 
changes across the kick - x remains constant because Px 
is absent from the coupling. That said, our approach is 
the following: 

(i) we first define the action-angle variables just before 
the n*'* kick: /„,6'„; 

(ii) action and angle are then transformed to the 
original cartesian variables via (a;„, Pxn) — 

{y/UZsinOn, \fnZcOs6n)\ 

(iii) changes across kick n are calculated in the cartesian 
coordinates via Aa; = 0, Apa; = kAsinikxn)', 

(iv) immediately after the kick we write Px^ = Pxn '^ 
Apx and x^ — Xn + Ax and change back to action- 
angle variables via /+ = (p^+ + a;+ )/2, 6'+ = 
arctan(x+/pa;,t)- 

(v) The last step is to propagate the preceding state to 
that instant n + 1 before the next kick via In+i — 



I, 



„+l, ^n+l 



e+ + T/./TT2i, 



n+l- 



Although the final composite map could be expressed 
with similar degree of complexity in terms of cartesian 
coordinates, use of action-angle variables is more conve- 
nient in view of the fact that action is conserved in the 
absence of perturbation. 

The final result can be cast in the form of an explicit 
map relating the dynamical states at kicks n and n+l: 



J 



1 -2/1 

In+1 = -Z I 2/„ sin 6rt 



-ek sm{k\/2I„ sin 0„) 



(3) 



9n+i — arctan 



2y/2I^cos6n + eksm{ky/2I^sm6n)J ^1 + 2I„+i 



T 



(4) 



where for convenience we introduce e = 2 A. The map has the noticeable feature of being fully explicit - from 



Eq. (jS]), In+i in Eq. ((4]) can be written in terms of 
/„ and 6n. In addition the map Jacobian has unitary 
norm, which guarantees the sympletic character of the 
theory. We also note that, in contrast to the standard 



J 



case, the map displays strong nonlinear dependence on 
the wave amplitude. A linearized map, obtainable only 
when e ^ 1, can be cast into the canonical form 



/„ = /„+i - {l/V2)ek^/I^cos9nCos{k ^/2T^s'me„) 
T ek sin 6'„ cos(fc ■\/2/„+i sin 0„) 



7n+l — 0,1 



^/lT2I, 



n+1 



ysZ 



(5) 
(6) 



n+l 



r 



The linear form can be naively obtained from the Hamil- 
tonian ([2]) alone if across the kick one assumes 9 — >■ 6'„ 
and / — > In+i'i our theory shows why these assumptions 
should be this way. 



III. ANALYTICAL AND NUMERICAL 
RESULTS 

We are now in position to explore the properties of the 
fully nonlinear map (|31) and ^. When e — >■ the action 
/ becomes constant and the phase 9 advances steadily 
according to 



T 



7n+l 



VTT2T 



(7) 



As mentioned earlier, we shall look into cases where both 
periods are similar, since a variety of cases fall in this 
ordering category. Let us further narrow our focus on 
the main resonance, the one for which each wave cycle 
corresponds to a full orbital gyration of the magnetized 
particles. This case corresponds to take 9n+i — 9n = 27r, 
which sets the position of the main resonance at 



/., 



87r2 



(8) 



Expression ([8]) is only an estimate that does not include 
any effects resulting from the perturbing wave. Never- 
theless, it provides a first step to investigate the role of 
resonances in the dynamical system. Recalling that the 
wave frequency is measured in units of the cyclotron fre- 
quency qBo/m, when wave and cyclotron frequencies co- 
incide, T — 27r/l = 2tt and Ires -^ 0. For larger values of 
the wave period T, the resonance moves towards higher 
values of the action, and for smaller values the resonance 
cannot be realized for positive values of the action. Im- 
portantly, when relativity is suppressed with / — >■ in 
Eq. d?]) , the phase advance becomes independent of the 
action and one falls in a degenerate theory where the res- 
onant island is absent Q . This feature contrasts with the 
relativistic standard case where suppression of relativistic 
effects does not remove the main resonance [ll|, [l2] ■ 



We now look into a case where the resonance is present 
in the phase-space with the period T slighly larger than 
the cyclotron period: T = 27r(l -I- 1/15) > 27r. We also 
consider fc = 2 as a representative wavevector of the 
modes analyzed, and display two distinct situations in 
Fig. [TJ ek^ = 0.1 in panel (a) and efc^ — 0.83 in panel 
(b). It will become clear later why we work with the 
combination efc^ instead of the wave amplitude e alone. 
In both cases the phase-spaces are fairly regular, but a 
remarkable fact distinguishes the panels apart. While 
in (a) the resonance takes the usual pendulum-like shape 
with elliptic and hyperbolic points approximately located 
at the same level along the action axis, in panel (b) one 
sees that the hyperbolic points move down to / ~ while 
the elliptic point stays at the same original location. The 
relevance of this feature lies in the possibility of parti- 
cle acceleration from low initial energies. Indeed, while 
in panel (a) particles launched with J « remain with 
small values of the action, in panel (b) one sees that par- 
ticles launched under the same condition / w perform 
much larger excursions along the island separatrix. In- 
terestingly, the linearized map dS]), ^ yields the angular 
position of the two hyperbolic points exactly at 7r/2 and 
3tt/2 for all values of e. Close examination of panel (b) 
shows however that this is only approximate, as one sees 
that both points sit to right of coordinates tt/2 and 3tt/2 
respectively. The region near 9 = 7r/2, in particular, is 
detailed in panel (c) where one can see the positioning 
of the hyperbolic point - numerical analysis reveals that 
it is located at 6* = 1.13 x 7r/2, / = 0.00157. The intrin- 
cated nonlinear dependence of the full map ([3]), (01) on 
£ displaces the hyperbolic points from their approximate 
phases, which is something to consider if one wishes to 
adjust wave-particle phase for optimum acceleration. 

Topology of Fig. [T]is similar in related contexts involv- 
ing time dependent but spatially uniform perturbations 
[20], although the fully analytical map dU, ^ reveals 
that the hyperbolic points in both cases do not coincide. 

Fig. [Ijb) also suggests that one way to improve the ac- 
celeration efficiency would be to increase the mismatch 
between cyclotron and wave frequency. According to the 
estimate of Eq. ([8]) , larger wave periods would pull the 
resonance upwards, and if one could choose a sufficiently 
large e, the hyperbolic points could still be brought down 




FIG. 1: Nonlinear resonance: pendulum-like with ek = 0.1 
in (a), and ek — 0.83 in (b). Panel (c) details the positioning 
of the hyperbolic fixed point of panel (b) near ■k/2. k — 2 and 
T = 2tt{1 + 1/15) in all cases. 



to J « enabling particles to coherently loop from small 
to large values of tlie action. There is however a natural 
limitation to this outlined procedure, which arises from 
the loss of stability of elliptic points as the wave ampli- 
tude e grows. In the case under study, one can expect 
the accelerating mechanism to work satisfactorily up to 



the point where the central point at tt looses stability. 
Beyond that, one can expect the intrusion of appreciable 
chaotic activity into the system with the concomitant loss 
of orbital coherence. It is true that even before the ellip- 
tic point bifurcates chaos is present, but for the sake of 
simplicity we take the bifurcation as indicating the limits 
of regular regimes in the system. The condition we seek 
for regular acceleration can be thus stated as the one 
associated with periods T for which one can maneuver 
the control parameters such that the hyperbolic points 
touch the axis 1 = before the elliptic point bifurcates. 
Maximum, or optimum acceleration can be achieved for 
the particular period T = Topt where touch down and 
bifurcation occur simultaneously. 

To obtain the parameters for which hyperbolic {h) 
points reach / = we look for fixed points of Eq. ([4]) 
with periodicity 27r at / = 0, i.e., On+i = On + 27r with 
6*,, = 6'„(mod27r) : 



T + arctan 



tan(6'ft) 



1 



efc2 tan(6'/i) 



= Oh + 27r. 



(9) 



As mentioned earlier, fixed points of Eq. ([9]) are not 
aligned with 9 = 7r/2, 37r/2 as in the linearized version of 
map ^ and Q. Nevertheless, careful inspection of Eq. 
([9]) along with the use of trigonometric identities allows 
to find the following solution for the fixed points: 



J 



tan 9h 



ek^ tan(T) ± ^tan(r)^e2fc4 tan(r) + Se fc2 - 16tan(T) 
2(£fc2-2tan(r)) ' 



(10) 



r 



Real solutions of expression ([T0|) exist only when 
e2fc4tan(T) + 8eP - 16tan(T) > 0, from which one ob- 
tains the corresponding threshold for touch down (td): 



{ek%d = 



tan(T) 



l + tan2(r)-l 



(11) 



When T = 27r(l -t- 1/15), as used in Fig. [T]for instance, 
{£k'^)td = 0.85066; this is why we chose that particular 
value of the product ek"^ in Fig. [T] (b) - it is close to the 
touch down limit. 

As mentioned, one has regular acceleration only when 
touch down occurs before the central fixed point at 9 = tt 
looses stability via a period doubling bifurcation (pd). 
To obtain the conditions for period doubling we first of 
all observe that analysis of Eq. Q at = tt allows to 
conclude that result ([5]) is in fact exact, independently 
of the value of e. Then, linear stability analysis of Eq. 
dl]) at the fixed point readily indicates that bifurcation 
occurs at: 



{£k'^)pd 



AT' 



TT {T'^ ~ At:'')' 



(12) 



We point out that from expressions (fTTj) and (|12|) it be- 
comes apparent that the effective parameter controlling 
the perturbative strength is the combination ek' and not 
the wave amplitude alone, as commented earlier. 

The next step is to compare both thresholds furnished 
by Eqs. dTT]) and dH]). This is done in Fig. [2] where 
one sees that there is indeed a critical point in the plane 
ek' versus T, where touch down and period doubling oc- 
cur simultaneously. For values of the period T below 
the one at the critical point, touch down occurs earlier 
than doubling, and after the critical point the ordering 
is reversed. At the critical point one attains optimum 
conditions for regular acceleration, as discussed earlier. 
The critical point period can be numerically evaluated 



as T, 



opt ~ 1.24 X 27r. We also have {ek')opt = 3.67843 
This means that the dimension- 

1.03 



and lopt = 0.26469. 



less momentum excursion reads Px ^ ^j2{2Iopt) 
which provides a significant amount of acceleration Vx = 
Px/ \J^ +Px ~ O-^S especially when one looks at ion dy- 
namics. To obtain the estimate for maximum momen- 
tum we simply double the value of the action at the fixed 



point, as suggested by the topology of Fig. [T] (b), and 
make use of expression (jS]) along with the relation be- 
tween action-angle and cartesian coordinates. Sopt ~ 1, 
so we are beyond the validity of linear approximation ([5]) , 



period— doubling 




er 



0.05 0.10 0.15 0.20 0.25 0.30 
T/27T - 1 

FIG. 2: Comparison of the threshold curves for period- 
doubling and touch down. Labels respectively indicated. The 
black dot indicates the point of optimum acceleration. 



IV. CONCLUSIONS 



wave amplitude, which affects the positioning of fixed 
points. Also, special regular orbits were shown to exist 
which can be used to accelerate particles. These special 
orbits bifurcate to chaos when the interaction parameter 
grows beyond a certain threshold. The orbits and all the 
associated bifurcation process are only present when non- 
linear relativistic mass correction of particles is properly 
taken into account. If relativistic effects are erroneously 
ignored, resonant islands are replaced with KAM curves 

Comparing ours with similar settings, we initially re- 
call that topology of Fig. [1] resembles that of related 
contexts invol ving time dependent but spatially uniform 
perturbations [2^ , although the fully analytical map ([3]) , 
Q reveals different positioning of fixed points. We point 
out as well that relativistic generalizations of the Kar- 
ney problem applied to the interaction of magnetized 
electrons and perpendicular electrostatic lower hybrid 
waves have also been studied, with closer focus either 
on pendulum-like islands or on waves with much larger 
than the relativistic cyclotron frequencies ;18]. Our prob- 
lem here is more directed to the resonant interaction and 
generation of the large distorted island. 



To conclude, we developed and investigated an analyt- 
ical canonical map describing the relevant case of impul- 
sive interaction of waves and magnetized relativistic par- 
ticles. In contrast to the better known standard canoni- 
cal case, the present map is shown to be nonlinear in the 
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